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The  Accuracy  of  Remapping  Irregularly  Spaced  Velocity  Data  onto  a  Regular  Grid 
and  the  Computation  of  Vorticity 

R.K.  Cohn  and  M.M.  Koochesfahani 


Abstract.  The  velocity  data  obtained  from  Molecular  Tagging  Velo- 
cimetry  (MTV)  are  typically  located  on  an  irregularly  spaced  meas¬ 
urement  grid.  To  take  advantage  of  many  standard  data  processing 
techniques,  the  MTV  data  need  to  be  remapped  onto  a  grid  with  a 
uniform  spacing.  In  this  work,  accuracy  and  noise  issues  related  to 
the  use  of  a  least-squares-fit  to  various  low  order  polynomials  for  the 
remapping  of  these  data  onto  a  uniformly  spaced  grid  and  the  subse¬ 
quent  computation  of  vorticity  from  these  data  are  examined.  This 
information  has  relevance  to  PIV  data  processing  as  well.  As  noted 
by  Spedding  and  Rignot  (1993),  the  best  estimate  of  the  velocity 
vector  acquired  through  the  use  of  tracer  techniques  such  as  PIV,  is  at 
the  midpoint  of  the  displacement  vector.  Thus,  unless  special  care  is 
taken,  PIV  data  are  also  initially  obtained  on  an  irregular  grid. 

As  in  the  results  of  Fouras  and  Soria  (1998),  the  error  in  the  re¬ 
mapped  velocity  and  the  calculated  vorticity  field  is  divided  into  a 
mean  bias  error  and  a  random  error.  In  the  majority  of  cases,  the 
mean  bias  error  is  a  more  significant  source  of  error  than  the  more 
often  quoted  random  error.  Results  of  the  simulation  show  that  the 
best  choice  for  remapping  is  the  use  of  a  least-squares  fit  to  a  2nd 
order  polynomial  and  the  best  choice  for  vorticity  calculation  is  to 
use  a  4th  order  accurate,  central,  finite  difference  applied  to  uniformly 
sampled  data.  The  actual  value  of  the  error  depends  upon  the  data 
density  and  the  radius  used  for  the  selection  of  velocity  measure¬ 
ments  to  be  included  in  the  remapping  process.  Increasing  the  data 
density  and  reducing  the  fit  radius  improve  the  accuracy. 

1  Introduction 

In  recent  years,  many  researchers  have  made  use  of  full-field, 
two-component  optical  velocity  measurement  techniques,  such  as 
Particle  Image  Velocimetry  (PIV),  to  derive  flow  quantities  such  as 
the  out-of-plane  vorticity  from  velocity  data.  The  velocity  field  ac¬ 
quired  from  PIV  is  normally  thought  to  be  gathered  on  a  uniformly 
spaced  grid  which  allows  for  a  variety  of  standard  post-processing 
methods  to  be  utilized.  The  development  of  Molecular  Tagging  Ve¬ 
locimetry  (MTV)  has  placed  an  additional  complication  on  the  cal¬ 
culation  of  flow  variables  in  that  the  data  are  not  normally  collected 
on  a  uniformly  spaced  grid.  This  paper  deals  with  the  questions  re¬ 
lated  to  remapping  MTV  data  onto  a  regularly  spaced  grid  and  the 
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methods  used  to  compute  the  out-of-plane  vorticity  component  from 
these  remapped  data  sets. 

Molecular  Tagging  Velocimetry  is  a  full-field  optical  diagnostic 
which  allows  for  the  non-intrusive  measurement  of  a  fluid  velocity 
field.  This  technique  takes  advantage  of  molecules  which  have  long- 
lived  excited  states  when  tagged  by  a  photon  source.  This  technique 
can  be  thought  of  as  the  molecular  equivalent  of  Particle  Image  Velo¬ 
cimetry.  Rather  than  tracking  particles  placed  in  the  flowing  me¬ 
dium,  the  luminescence  of  regions  of  the  flow  containing  the  tracer 
molecules  is  tracked.  A  more  complete  description  of  the  imple¬ 
mentation  of  the  molecular  tagging  technique,  its  applications,  and 
the  parameters  necessary  for  an  optimal  experiment  can  be  found  in 
Gendrich  and  Koochesfahani  (1996)  and  Gendrich  et  al.  (1997),  and 
Koochesfahani  (1999).  The  accuracy  of  velocity  measurements  made 
using  Molecular  Tagging  Velocimetry  is  comparable  to  the  digital 
version  of  PIV  (DPIV). 

In  the  implementation  of  MTV,  a  series  of  laser-lines  is  used  to 
generate  a  two-dimensional  spatial  distribution  in  the  intensity  field 
within  the  flowing  medium.  Velocity  vectors  are  calculated  at  the 
intersection  of  these  laser-lines.  Generally,  the  measurement  loca¬ 
tions  are  not  uniformly  spaced.  Thus,  it  is  necessary  to  place  the 
velocity  data  onto  a  regular  grid  before  flow  variables,  such  as  vor¬ 
ticity,  can  be  computed  via  standard  finite  difference  techniques.  It 
should  be  noted  that  even  though  it  is  possible  to  generate  a  series  of 
regularly  spaced  laser-lines  in  the  flow,  it  is  still  necessary  to  remap 
the  data.  Both  an  unpublished  study  conducted  at  Michigan  State 
University  and  Spedding  and  Rignot  (1993)  have  reported  that  the 
best  estimate  of  the  location  of  the  velocity  vector  determined  by  a 
measurement  technique  which  tracks  a  tracer  in  a  flow  is  located  at 
the  midpoint  of  the  displacement  vector.  Thus,  unless  special  care  is 
taken  in  the  selection  of  the  measurement  windows,  data  collected 
from  PIV  measurements  is  not  on  a  uniformly  spaced  grid. 

Few  studies  have  examined  the  effect  of  remapping  randomly 
spaced  velocity  data  onto  a  regular  grid.  Agui  and  Jimenez  (1987) 
reported  that  low  order  polynomial  fits  and  “kriging”  techniques 
produced  the  most  accurate  representation  of  the  actual  velocity  field. 
However,  the  advantage  was  small  with  respect  to  other  methods  and 
no  quantitative  information  on  the  performance  of  the  polynomial 
and  kriging  methods  was  given.  Spedding  and  Rignot  (1993)  com¬ 
pared  an  inverse  distance  approach  with  the  use  of  a  “global  basis 
function”  and  found  that  the  global  basis  function  produced  generally 
more  accurate  results;  however,  the  results  were  highly  dependent 
upon  the  measurement  density. 

Several  authors  have  examined  the  accuracy  of  various  means  to 
compute  vorticity  from  velocity  data  already  on  a  regular  grid.  Sped¬ 
ding  and  Rignot  (1993)  used  a  1st  order  finite  difference  technique  for 
the  inverse  distance  method  and  directly  differentiated  the  global 
basis  function  to  compute  vorticity.  It  was  found  that  direct  differen¬ 
tiation  of  the  global  basis  function  produced  generally  superior  re¬ 
sults.  However,  as  with  the  velocity  results,  the  accuracy  was  highly 
dependent  upon  the  ratio  of  a  characteristic  length  scale  of  the  flow, 
L,  to  the  mean  spacing  between  measurements,  5.  Abrahamson  and 
Lonnes  (1995)  found  that  calculating  vorticity  by  computing  the  local 
circulation  around  a  point  resulted  in  slightly  more  accurate  vorticity 


results  than  differentiating  a  least-squares  fit  to  a  model  velocity 
field.  Luff  et  al.  (1999)  compared  the  1st  and  2nd  order  finite  differ¬ 
ence  methods  and  an  8-point  circulation  method  in  the  calculation  of 
vorticity  in  the  presence  of  both  noise  and  missing  data  points.  In 
terms  of  only  the  computed  vorticity  RMS,  the  1st  order  finite  differ¬ 
ence  technique  produced  the  best  results. 

One  shortcoming  of  the  above  mentioned  studies  is  that  only  the 
random  component  of  the  error  field  is  examined.  Fouras  and  Soria 
(1998)  found  that  the  error  in  the  vorticity  field  could  be  better  repre¬ 
sented  if  it  is  divided  into  two  portions:  a  mean  bias  error  due  to  spa¬ 
tial  filtering,  and  a  random  error  resulting  from  the  propagation  of 
error  in  the  velocity  measurements  into  the  vorticity  calculation.  In 
some  cases,  the  mean  bias  error  can  be  significantly  larger  than  the 
random  error.  This  study  recommends  differentiating  a  2nd  order 
polynomial  least-squares  fit  to  the  velocity  data  for  the  calculation  of 
vorticity  based  on  the  21  closest  points.  However,  at  low  data  densi¬ 
ties,  this  produces  larger  bias  errors  then  the  use  of  a  finite  difference 
method.  The  results  based  on  differentiating  the  fit  were  sensitive  to 
the  number  of  points  used  in  the  fit.  This  work  is  based  entirely  on 
regularly  sampled  velocity  data;  issues  connected  to  remapping  an 
irregular  data  set  were  not  considered. 

The  aforementioned  investigations  suggest  different  optimum 
methods  for  vorticity  computation  depending  upon  the  criterion  used 
to  assess  the  error.  In  our  work,  we  directly  compare  several  of  the 
different  vorticity  calculation  methods  which  were  determined  in  the 
previous  studies  to  produce  the  best  results.  In  addition,  we  include 
the  effect  of  the  remapping  of  the  velocity  field  on  the  estimation  of 
the  vorticity  is  also  considered. 


2  Comparison  Method 


The  present  study  makes  use  of  a  simulation  of  an  Oseen  vortex 
in  order  to  study  the  effect  of  remapping  an  irregularly  sampled  ve¬ 
locity  field  onto  a  regular  grid  and  the  calculation  of  the  out-of-plane 
vorticity  component.  This  flowfield  has  also  been  used  in  the  works 
of  Spedding  and  Rignot  (1993),  Fouras  and  Soria  (1998),  and  Luff  et 
al.  (1999).  The  azimuthal  velocity,  ue,  and  out-of-plane  vorticity,  o>z 
of  this  flow  field  are  described  by: 
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In  order  to  simulate  the  irregular  sampling  found  in  the  original 
velocity  field  measurements,  the  simulation  data  are  irregularly 
spaced.  The  irregular  spacing  is  generated  by  sub-dividing  the  meas¬ 
urement  field  into  8x5  sized  regions,  where  8  is  the  mean  spacing 
between  velocity  measurement  points.  A  random  number  generator 
is  then  used  to  determine  a  location  for  the  simulated  velocity  within 
each  5x8  sized  region.  Equation  (1)  is  then  used  to  establish  the 
velocity  at  this  location.  In  this  manner,  the  mean  spacing  between 
measurement  points  remains  equal  to  5;  however  the  actual  location 
of  the  measurement  varies. 

The  random  error  inherent  in  MTV  and  PIV  measurements  is 
simulated  by  the  addition  of  noise  to  the  velocity  field.  The  method 
used  is  similar  to  that  in  Luff,  et  al.  (1999).  A  random  number  gen¬ 
erator  is  used  to  add  a  random  percentage  of  noise,  with  a  maximum 
value  of  n%  to  each  component  of  the  cartesian  components,  i.e.  u 
and  v,  of  the  velocity  field.  Using  this  formulation,  the  velocity  at 
each  point  in  the  simulation  has  a  value  of: 

U  ~  U act  (1  +  n  random  )  (2) 

v  =  nc,( 


where  nrandom  is  a  random  number  with  a  value  -n  <  nrandom  <  +n . 


The  quantities  uact  and  vact  represent  the  actual  cartesian  velocities 
determined  from  we  in  equation  (1).  Although  error  values  ranging 
from  0%  <  n  <  10%  were  examined  in  the  simulations,  only  0%  and 
6%  values  are  presented  as  they  are  representative  of  the  other  noise 
values. 

The  velocity  data  in  the  present  study  are  remapped  onto  a 
regular  grid  by  means  of  a  local  least-squares  fit  to  a  two-dimensional 
2nd,  3rd,  or  4th  order  polynomial.  The  u  and  v  velocity  fields  are  fit 
separately.  Only  the  velocity  measurements  located  within  the  fit 
radius,  i?,  from  the  regular  grid  point  are  used  in  the  fitting  procedure. 
In  this  study,  this  radius  is  normalized  by  the  mean  spacing  between 
velocity  measurements,  8.  For  all  cases,  the  number  of  points  used  in 
the  fit  is  such  that  the  least-squares  fit  is  over-determined.  After  the 
fits  for  the  two  velocity  components  are  generated,  each  of  the  local 
fits  is  evaluated  at  the  coordinate  of  the  regular  grid  point  in  order  to 
determine  the  velocity  at  that  location.  Clearly  choosing  a  value  of  R 
that  is  too  large  will  result  in  a  large  amount  of  spatial  filtering  of  the 
data,  while  a  small  value  of  R  will  limit  the  ability  of  the  fit  to  reduce 
the  random  noise  present  in  the  original  data.  Note  that  the  order  of 
the  polynomial  places  a  limit  on  the  minimum  size  of  R  that  can  be 
used.  For  all  of  the  studies  conducted,  the  density  of  the  remapped, 
uniformly  spaced  grid  remains  the  same  as  the  initial  irregularly 
spaced  measurement  grid. 

Four  methods  are  used  to  estimate  the  out-of-plane  vorticity,  i.e, 


(3) 
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The  first  two  methods  estimate  the  two  derivatives  in  the  definition  of 
vorticity  by  means  of  a  1st  or  2nd  order  central  finite  difference  tech¬ 
niques  (2nd  and  4th  order  accurate,  respectively).  The  third  method 
performs  a  direct  differentiation  of  the  polynomial  least-squares  fit 
used  in  the  remapping  of  the  velocity  field.  This  method  has  the 
advantage  that  it  can  be  used  to  estimate  the  vorticity  at  any  point 
within  measurement  region.  The  final  vorticity  calculation  method 
computes  the  circulation  of  the  8  points  in  the  rectangular  region 
extending  one  regular  grid  point  in  each  direction  around  the  point  to 
be  examined.  The  calculated  circulation  value  is  then  divided  by  the 
area  in  order  to  determine  the  vorticity.  This  method  has  been  shown 
in  Raffel  et  al.  (1998)  to  be  identical  to  a  filtered  version  of  the  1st 
order  central  finite  difference  technique.  Figure  1  illustrates  the  data 
points  used  for  the  calculation  of  the  vorticity  in  these  various  meth¬ 
ods. 


This  study  examines  the  effect  of  varying  the  mean  data  density, 
L/ 5  and  the  normalized  fit  radius,  RJ 5  on  the  accuracy  of  the  remap¬ 
ped  velocity  field  and  the  calculated  vorticity  field.  The  characteris¬ 
tic  flow  scale,  L,  used  in  this  study  is  the  vortex  core  radius,  rcore 
defined  as  the  distance  from  the  peak  vorticity  to  the  location  where 
the  vorticity  has  dropped  by  a  factor  of  e1 .  Simulations  are  con¬ 
ducted  for  values  of  L/ 5  ranging  from  2  to  10  and  for  Rj 8  ranging 
from  2  to  6.  Note  that  in  order  for  the  fit  to  be  over-determined,  not 
all  of  these  values  can  be  used  for  all  orders  of  the  least-squares  fit. 

As  in  the  results  of  Fouras  and  Soria  (1998),  the  accuracy  of 
both  the  velocity  and  vorticity  calculations  methods  are  assessed  in 
terms  of  the  mean  bias  error  caused  by  spatial  filtering  and  the  ran¬ 
dom  error.  Both  the  propagation  of  the  error  in  the  original  meas¬ 
urements  to  the  remapped  field  and  the  placement  of  the  randomly 
spaced  points  onto  the  regular  grid  generate  the  random  error.  Within 
each  sample  set,  100  independent  simulations  are  conducted.  The 
mean  bias  error  will  be  denoted  by  the  subscript  bias  and  refers  to 
difference  between  the  mean  value  of  these  100  velocity  (or  vorticity) 
measurements  and  the  exact  value  at  each  point  in  the  flowfield  de¬ 
termined  from  equation  (1).  The  random  error  is  quantified  by  the 
rms  of  the  velocity  (or  computed  vorticity)  data  in  the  sample  set. 
For  example,  for  the  x-component  of  velocity,  these  are  defined  as: 


i  100 


1  100 

u = — y  u. . 
100“  ' 


(4) 


In  these  expressions,  u{  is  the  velocity  at  a  particular  point  of  an  indi¬ 
vidual  realization.  All  of  the  velocity  and  vorticity  error  values  re¬ 
ported  here  are  normalized  by  their  respective  maximum  values  de¬ 
termined  from  equation  (1).  The  normalized  values  are  shown  with¬ 
out  the  ( )*. 

3  Remapping  Results 

Figure  2a  displays  the  mean  velocity  bias  error  for  the  2nd,  3  rd, 
and  4th  order  polynomial  fits.  The  results  presented  in  this  figure  are 
for  the  case  of  0%  added  noise  because  the  addition  of  noise  has  no 
effect  on  the  bias  error,  which  is  a  mean  quantity.  For  all  three  fits, 
the  mean  velocity  data  density  is  kept  fixed  at  Lib  =  3.5.  That  is, 
there  are  nominally  7  velocity  vectors  along  the  vortex  core  diameter. 
For  the  2nd  and  3rd  order  polynomial  fits,  three  different  values  of  RJ 5 
are  examined.  However,  only  the  two  larger  values  are  used  for  the 
4th  order  polynomial  fit  to  ensure  there  are  enough  data  points  avail¬ 
able  for  the  fit.  In  terms  of  the  mean  bias  error,  reducing  this  radius 
results  in  a  significant  decrease  in  the  bias  error.  For  the  2nd  and  3rd 
order  polynomial  fits,  reducing  from  RJ 5  =  4  to  RJ 5  =  2  results  in  a 
decrease  of  the  peak  mean  bias  error  from  8%  to  less  than  1%.  This 
effect  is  present,  although  less  dramatic  in  the  results  for  the  4th  order 
polynomial  fit. 

In  terms  of  the  mean  bias  error,  the  most  accurate  results  are 
obtained  using  the  least-squares  fit  to  a  4th  order  polynomial.  How¬ 
ever,  the  difference  in  the  bias  error  between  the  4th  order  fit  using 
RJ 8  =  3,  and  the  2nd  order  fit  with  RJ 8  =  2  is  approximately  0.4%.  It 
is  interesting  to  note  that  the  results  for  the  2nd  and  3rd  order  polyno¬ 
mial  fits  are  nearly  identical.  For  all  three  fit  orders,  the  peak  bias 
error  occurs  at  approximately  0.6rcore .  In  the  region  r/rcore  >  1.5 ,  the 
velocity  values  tend  to  be  overestimated,  rather  than  underestimated. 
Since  the  area  of  the  region  r/rcore  >  7.5  is  significantly  larger  than 
the  region  where  the  velocity  values  are  underestimated,  the  net  aver¬ 
age  of  the  bias  error  becomes  very  small.  Thus,  one  should  be  cau¬ 
tious  about  the  use  of  an  accuracy  measure  which  is  averaged  over 
the  entire  structure  as  this  does  not  represent  the  actual  error  seen  at 
any  individual  measurement  location. 

Figures  2b  and  2c  show  the  random  error  found  in  the  velocity 
remapping  field  for  cases  of  0%  and  6%  added  noise  respectively. 
From  Figures  2b  and  2c  it  can  be  seen  that,  generally,  the  value  of  the 
rms  error  is  less  than  2%  at  all  locations.  This  value  can  only  be 
reduced  by  a  maximum  of  1.5%  by  the  optimal  choice  of  fit  order  and 
R/b,  whereas  a  reduction  of  8%  is  seen  in  the  bias  error.  It  is  also 
interesting  to  note  that  in  the  0%  added  noise  results,  reducing  R/b 
results  in  a  decrease  in  the  random  error.  However,  for  the  6%  added 
noise  results,  reducing  R/b  results  in  an  increase  in  the  random  error. 

Unless  otherwise  stated,  the  remainder  of  the  results  presented  in 
this  paper  will  use  the  2nd  order  polynomial  for  the  remapping  proc¬ 
ess.  This  choice  is  made  because  of  this  condition  has  a  very  small 
bias  error,  and  it  is  not  as  computationally  intensive  as  the  3rd  and  4th 
order  polynomials.  Figure  3  shows  the  effect  of  grid  density  on  the 
accuracy  of  the  remapping  for  the  2nd  order  polynomial  fit.  As  ex¬ 
pected  shows  that  increasing  the  grid  density  can  significantly  reduce 
the  bias  error  in  the  remapping.  In  order  to  achieve  a  bias  error  of 
less  than  1%,  LJ 5  should  be  greater  than  3.  Furthermore,  increasing 


the  grid  density  also  results  in  a  significant  decrease  in  the  random 
error. 

4  Vorticity  Calculation  Results 

In  this  section,  the  error  generated  by  the  four  methods  for  cal¬ 
culating  the  out-of-plane  vorticity  field  will  be  examined.  First,  we 
will  discuss  the  results  from  directly  differentiating  the  various  poly¬ 
nomial  orders  used  in  the  remapping  procedure  in  order  to  determine 
the  vorticity  value.  Then,  the  results  from  this  method  will  be  com¬ 
pared  with  those  from  the  finite  difference  and  circulation  methods 
applied  to  the  remapped  data  using  the  2nd  order  polynomial. 

Figure  4a  shows  the  mean  bias  error  for  differentiating  the  2nd, 
3rd,  and  4th  order  polynomials  for  several  values  of  R/b.  Similar  to  the 
velocity  bias  error  results,  decreasing  R/b  decreases  the  bias  error. 
The  smallest  values  for  the  mean  bias  error  are  found  using  the 
smallest  values  of  R/b  and  the  3rd  and  4th  order  polynomials.  These 
selections  result  in  a  bias  error  of  less  than  4%.  It  is  interesting  to 
note  that  although  the  error  in  the  remapping  of  the  velocity  field 
through  the  use  of  the  2nd  and  3rd  order  polynomials  are  nearly  identi¬ 
cal,  the  vorticity  estimates  by  differentiating  these  polynomials  differ. 

Figures  4b  and  4c  show  the  effect  of  the  order  of  the  polynomial 
fit  on  the  random  component  of  the  error.  For  the  case  in  which  no 
noise  is  added  to  the  velocity  data,  the  random  error  is  less  than  1.5%. 
The  addition  of  6%o  random  noise  results  in  a  small  increase  in  the 
random  component  of  the  error  for  the  majority  of  cases.  For  the  3rd 
order  fit  and  R/b  - 2, \  a  large  spike  is  seen  in  the  error.  It  is  believed 
that  this  spike  is  generated  because  the  number  of  points  utilized  for 
the  fit  is  only  slightly  larger  than  the  minimum  number  of  points 
required.  With  the  exception  of  that  case,  there  is  generally  little 
difference  in  the  random  error  among  the  various  polynomial  fits. 
Since  differentiating  the  4th  order  polynomial  with  7?/5  -  Jhas  one  of 
the  smallest  bias  errors  and  its  random  error  is  less  than  that  of  the  3rd 
order  fit,  only  the  vorticity  estimated  from  the  4th  order  fit  will  be 
compared  with  those  calculated  using  the  finite  difference  and  circu¬ 
lation  methods. 

Figure  5a  compares  ®bias  found  by  the  four  different  methods 
considered  here.  Once  again,  the  effect  of  adding  noise  to  the  initial 
velocity  field  on  the  bias  error  is  negligible,  therefore,  only  the  case 
of  0%  added  noise  is  shown.  The  qualitative  features  of  the  four 
methods  are  very  similar.  The  maximum  (nbias  occurs  at  r/rcore  -  0 
which  is  the  location  of  the  peak  vorticity.  For  r/rcore  >  1.5,  there  is  a 
small  overshoot  where  the  vorticity  value  is  overestimated.  Note  that 
although  the  numerical  amount  of  the  overshoot  is  small  relative  to 
that  of  the  undershoot,  the  area  occupied  by  the  region  of  overshoot  is 
roughly  three  times  larger  than  the  region  of  undershoot.  Thus,  the 
overall  area-averaged  vorticity  bias  error  is  very  small.  As  a  result, 
the  estimate  of  the  overall  circulation  of  the  vortex  computed  by  inte¬ 
grating  the  vorticity  field  from  any  of  these  methods  is  accurate  to 
better  than  0.1%  even  though  the  peak  bias  error  for  these  cases  is 
nearly  20%. 

In  terms  of  the  mean  bias  error,  it  is  apparent  that  differentiating 
the  4th  order  polynomial  fit,  and  the  use  of  the  2nd  order  finite  differ¬ 
ence  technique  produce  the  most  accurate  vorticity  field  information. 
For  comparison  purposes,  the  results  from  R/b  =  4  are  also  shown. 
As  expected,  reducing  the  value  of  R  tends  to  improve  the  accuracy 
of  the  vorticity  calculation  as  well  as  decreasing  the  difference  be¬ 
tween  the  accuracy  of  the  two  techniques.  The  circulation  method 
and  the  1st  order  finite  difference  method  produce  results  with  a  sig¬ 
nificantly  larger  bias  error  than  the  other  two  methods. 

Figure  5b  and  5c  show  the  random  error  for  the  vorticity  calcu¬ 
lation  methods  for  the  cases  of  0%  and  6%  added  noise  respectively. 
Generally,  differentiating  the  4th  order  polynomial  produces  the 


smallest  random  error  while  the  2nd  order  finite  difference  method 
produces  the  largest.  As  with  the  velocity  field,  the  improvement 
which  can  be  realized  through  the  use  of  the  optimal  method  to 
minimize  the  random  error  is  much  smaller  than  that  which  can  be 
realized  by  minimizing  the  bias  error.  However,  the  difference  in  the 
random  error  between  these  two  techniques  is  less  than  1%. 

As  expected,  increasing  the  density  of  the  original  data  also 
dramatically  reduces  the  mean  bias  error  as  seen  in  Figure  6a.  This 
figure  only  shows  results  for  the  2nd  order  finite  difference  technique 
and  RJ 5  =  2,  however,  the  qualitative  features  of  all  of  the  methods 
are  identical.  Increasing  the  mean  data  density,  L/b,  from  2.5  to  4 
results  in  a  decrease  of  the  mean  bias  error  from  20%  to  approxi¬ 
mately  1%.  In  terms  of  the  random  error,  increasing  the  density  of 
the  data  also  results  in  a  decrease  of  the  random  error. 

5  Conclusions 

The  effect  of  remapping  irregularly  spaced  velocity  measure¬ 
ments  onto  a  uniformly  spaced  grid  and  the  accuracy  of  the  out-of- 
plane  vorticity  computed  from  this  information  are  studied  through 
the  use  of  a  Gaussian  core  vortex  simulation.  Remapping  onto  a 
regular  grid  was  performed  by  the  use  of  a  least-squares  fit  to  2nd,  3rd, 
and  4th  order  polynomials.  Four  methods  are  used  for  the  calculation 
of  vorticity  from  the  velocity  data.  The  derivatives  necessary  for  the 
computation  of  vorticity  can  be  found  by  directly  differentiating  the 
polynomials  determined  from  the  least-squares  fit  or  by  performing  a 
1st  or  2nd  order  finite  difference  calculation  on  the  regularly  spaced 
data.  The  final  method  computes  the  local  circulation  of  the  region 
around  a  point  and  divides  this  circulation  by  the  area.  The  effect  of 
varying  the  normalized  grid  density  (ratio  of  the  flow  characteristic 
length  to  the  mean  spacing  in  the  initial  velocity  measurement)  and 
the  maximum  radius  from  which  points  are  used  in  the  remapping 
process  are  examined.  For  all  of  the  studies  conducted,  the  density  of 
the  remapped,  uniformly  spaced  grid  remains  the  same  as  the  initial 
irregularly  spaced  measurement  grid. 

As  in  the  study  by  Fouras  and  Soria  (1998),  the  error  resulting 
from  the  remapping  and  the  calculation  of  the  out-of-plane  vorticity 
is  divided  into  a  mean  bias  error  due  to  spatial  filtering  and  a  random 
error  due  to  the  remapping  process  itself  and  the  propagation  of  the 
error  in  the  original  data.  Generally,  the  errors  resulting  from  the  bias 
error  are  significantly  larger  than  random  error.  The  mean  bias  error 
is  not  affected  by  noise  in  the  original  measurements  and  can  be  de¬ 
creased  by  increasing  the  grid  density.  The  random  error  is  affected 
by  the  presence  of  noise  in  the  original  measurements,  which  causes 
an  increase  in  the  random  error.  However,  the  filtering  inherent  in 
the  least-squares  fitting  process  tends  to  decrease  the  magnitude  of 
the  random  error  in  the  remapped  data. 

In  terms  of  the  mean  velocity  bias  error,  it  is  necessary  for  the 
grid  density  to  be  suitably  high  and  for  the  fit  to  be  local,  i.e.,  small 
values  of  RJ 5,  to  generate  an  accurate  remapping.  In  this  study,  the 
4th  order  polynomial  produced  the  most  accurate  remapping;  how¬ 
ever,  due  to  the  large  number  of  points  needed  for  the  fit  to  be  deter¬ 
mined,  it  is  not  possible  to  use  the  same  small  values  of  R/b  as  can  be 
used  with  the  lower  order  polynomials.  The  difference  in  the  bias 
error  between  the  4th  order  with  RJ 5  =  3  and  the  2nd  order  with  R/b  = 
2  is  very  small.  Thus,  it  is  felt  that  the  use  of  a  2nd  order  polynomial 


for  the  remapping  is  appropriate.  In  order  to  obtain  results  with  a 
maximum  mean  velocity  bias  error  of  approximately  1%  or  less  for 
the  2nd  order  polynomial,  the  normalized  data  density  L/b  should  be 
greater  than  3,  and  R/b  should  be  less  than  or  equal  to  3.  The  differ¬ 
ence  in  the  random  component  of  the  error  between  the  various  fit 
orders  is  relatively  small  compared  to  the  changes  seen  in  the  bias 
error.  The  rms  among  the  methods  only  varied  by  approximately  2% 
of  the  peak  value  between  the  various  methods. 

The  most  accurate  results  for  the  vorticity  calculation  were  ob¬ 
tained  by  differentiating  the  4th  order  polynomial  fit  and  by  use  of  the 
2nd  order  finite  difference  method.  The  vorticity  results  are  more 
sensitive  to  the  choice  of  R/b  and  L/b  than  the  velocity  remapping 

results.  However,  L/b  >  3.5  and  R/b  <2.5  yields  vorticity  values 
with  both  a  mean  bias  and  random  error  of  less  than  approximately 
2%.  It  should  be  noted  that  the  value  quoted  is  the  maximum  error 
across  the  vortex  core.  The  error  at  most  points  is  considerably  less 
than  2%.  As  with  the  velocity  error,  little  change  is  seen  among  the 
various  methods  in  terms  of  the  random  error. 

Although  the  accuracies  of  direct  differentiation  of  the  4th  order 
polynomial  and  the  use  of  the  2nd  order  finite  difference  method  are 
comparable,  it  is  believed  that  the  2nd  order  finite  difference  tech¬ 
nique  is  better  suited  for  the  calculation  of  vorticity  values  from  the 
PIV  or  MTV  data  since  this  method  is  less  computationally  intensive. 
It  should  be  noted  that  if  the  actual  value  of  R/b  and  L/b  are  known  in 
an  experiment,  it  is  possible  to  use  the  results  of  these  simulations  to 
extrapolate  and  estimate  the  actual  vorticity  value  from  under¬ 
resolved  data. 
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Figure  1 :  Velocity  measurement  locations  used  in  the  estimation  of  the  spatial  derivatives,  (a)  1st  order  finite  difference,  (b)  2nd  order  finite 
difference,  (c)  8-pt  circulation  method. 


Figure  2:  Accuracy  of  remapped  velocity  field  for  three  different  values  of  R/5.  (a)  Mean  bias  error,  (b)  Random  error  with  0%  noise 
added  to  initial  velocity  field,  (c)  Random  error  with  6%  noise  added  to  initial  velocity  field. 


Figure  3:  Accuracy  of  remapped  velocity  field  for  different  original  velocity  data  densities  for  the  2nd  order  polynomial  fit.  (a)  Mean  bias 
error,  (b)  Random  error  with  0%  noise  added  to  initial  velocity  field. 
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Figure  4:  Accuracy  of  out-of-plane  vorticity  field  computed  by  differentiating  the  local  polynomial  fit  for  different  values  of  R/3.  (a)  Mean 
bias  error,  (b)  Random  error  with  0%  noise  added  to  initial  velocity  field,  (c)  Random  error  with  6%  noise  added  to  initial  velocity  field. 


Figure  5:  Accuracy  of  out-of-plane  vorticity  field  computed  by  four  different  calculation  methods  for  different  values  of  R/S.  (a)  Mean  bias 
error,  (b)  Random  error  with  0%  noise  added  to  initial  velocity  field,  (c)  Random  error  with  6%  noise  added  to  initial  velocity  field. 


Figure  6:  Accuracy  of  out-of-plane  vorticity  field  for  different  original  velocity  data  densities.  Vorticity  is  calculated  using  the  2nd  order 
finite  difference  method,  (a)  Mean  bias  error,  (b)  Random  error  with  0%  noise  added  to  initial  velocity  field. 


